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1. Introduction 

Rigid body motion has always been an interesting and very usable subject of classical 
mechanics PUI21- Many subtle points of dynamics and of mathematical techniques for 
studying the behaviour of physical systems can be learned by studying the motion of 
rigid bodies. On this matter see, for example, [SI IU El El 13 El El EH] • Part of the interest 
of today comes from the insight that can be gained on the behaviour of spinning asteroids 
or artificial satellites and, furthermore, rigid body motion can be chaotic El ITH] . 
Besides, many everyday phenomena can be understood in terms of rigid body motion 
at least in a first approximation. For example, balls rolling on inclines, motion in 
toboggans, bowls in a bowling alley, the motion of snow boarders, the dynamics of 
bicycles or wheels, the behaviour of billiard balls, and so on. Though the importance 
of rigid bodies is clear, some of the problems involving rolling particles on a surface are 
often modelled in beginning courses as point particles sliding on surfaces [HI E] • This 
modelling is an appropriate pedagogical device in introductory courses but we want to 
show here that the problems can be addressed using a rigid body approach in more 
advanced courses. 

In this work we analyse the motion of a sphere that rolls without slipping on the 
inside of a right circular cone under the influence of a uniform gravitational field acting 
verically downwards, in the direction of the symmetry axis of the cone. The motion 
of a spherical body rolling without slipping on surfaces of revolution has been recently 
studied with the purpose of illuminating control processes ^0]- Here our aim is to 
study the motion of a sphere on the inner surface of a conical surfaceas an exactly 
solvable example of rigid body motion. We obtain the general solution of this problem 
expressing it in quadratures. We analyse certain qualitative features of the motion, like 
the existence of a stable circular orbit, establishing an analogy with particle motion in a 
central force field in two dimensions. For the case of a conical surface in free fall we find 
that the general solution of the problem can be casted in terms of expressions similar 
to those defining sets of ellipses and hyperbolas. We also calculate the apsidal (apogee) 
angle of the center of mass (CM) motion to find the orbit's symmetry axes and argue 
that the sphere's CM trajectory densely fills a strip on the conical surface. 

2. The equations of motion 

To describe the position of the sphere's centre of mass (CM), we choose a Cartesian 
coordinate system (x, y, z) such that the origin of coordinates is at the position of the 
sphere's centre of mass (CM) when the sphere simply rests on the cone (that is, it 
corresponds to the vertex of the imaginary cone — shown dashed in figure 1 — on which 
the CM moves); the coordinates are then (see Figure 1) 

x = r sin a cos??, 

y = r sin a sin (1) 

z = r cos a, 
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Figure 1. The cone on which the sphere of radius a rolls is shown in a continuous 
black line. The dashed cone is the imaginary cone on which the sphere's centre of mass 
(CM) moves. The imaginary cone is the same as the actual one but displaced upwards 
a distance a esc a. The origin of coordinates is chosen as the vertex of the imaginary 
cone so the vertex of the actual cone is at zq — —a esc a. The figure also shows the 
CM generalised coordinates r and The radius of the parallel circle containing the 
CM is symbolised by p. 



where a is half the angle of aperture of the cone — therefore, < a < tt/2, a = n/2 
corresponds to motion on a plane, whereas a = to motion on a cylindrical surface — , 
r is the distance from the origin to the CM, and t? is the polar angle of the sphere's CM. 
All these relations may be seen from figure 1. The tangential components of the CM 
velocity, in the, respectively, meridional and parallel directions to the conical surface 
(using a sort of geographical terminology), are 



where M is the mass, I = 5Ma 2 /2, is the moment of inertia of the sphere, a its radius; 
g is the acceleration of gravity, and <f), 9 and tp are Euler's angles pQ. Lagrangian (jH)) 
with u, v, and z defined appropriately, is fairly general; the description it gives is valid 
for an arbitrary surface of revolution with a sphere rolling on its surface. 

The motion of the sphere is further constrained according to u — auj ■ e\ = 0, 
and v + au ■ e 2 = 0, where uj is the angular velocity vector of the sphere, and 
ei = (sin a cos sin a sin cos a) , e 2 = (— sin??, cosi9, 0) are unit vectors in the, 
respectively, direction of the meridional line and the parallel circle on the cone. The 



u = r, 



(2) 



v = r sina$. 



Therefore, Lagrangian of the system can be written as 
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vector 10 can be written in terms of the Euler angles, an expression that can be found in 
many textbooks P3I2JIH]- If we employ equations (J2J) for substituting u and v, we obtain 



6 sin( 

a 



$) +ip sin 6 cos(0 — ■&) 



0. 



$ sin a — h 9 sin a cos(</> — $) — <f> sin a + 
a 

ip [sin a sin # sin (0 — i?) + cos a cos #] =0; (4) 

these two relations just mean that the sphere rolls with no slipping on the surface of 
the cone. 

The equations of motion of the sphere can be obtained, from L and constrictions 
(HJ), using Lagrange's undetermined coefficient method^], to get 



ii + v 



v — u 



sin a 



-clklu 



7 



g cos a, 



sin a 



au + uvk = 0, 



(5) 



where p is the radius of curvature of the parallel circle at point of contact: p = r sin a, 
K = cosa/p, and uj = is the component of u normal to the cone. 

The first two equations above are essentially Newton's second law for the centre of 
mass of the sphere, the third one is its rotating counterpart for a rotating rigid sphere. 
With the help of the constrictions, the Euler's angles have been eliminated from the 
description. 



3. Conserved Quantities in the Motion of the Sphere 

It should be clear that the total energy, E, of the sphere is a constant of motion since we 
are neglecting any dissipative effects. The energy of the sphere can be obtained directly 
from equations (jHJ), 

E = ^M(u 2 + v 2 ) + ^Ico 2 + 7 -Mgz. (6) 

We note in passing that this expression © is valid, again with u and v properly defined, 
for any surface of revolution. The point z = has been chosen as the point of zero 
potential energy PUEI]- It can be worth noting that the factor 7/10 in the CM's kinetic 
energy, instead of the usual 1/2, comes from effects associated with the rolling-but-no- 
slipping condition. 

For arbitrary surfaces of revolution, the energy is the only constant of motion. But 
in our conical case it is simple to show, using the equations of motion (jHJ), that the 
z-component of the angular momentum of the center of mass (CM) is also a constant 
of the motion, 

L z = Mpv = Mp 2 §. (7) 
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Moreover, for spheres rolling on flat surfaces — with symmetry axis as z-axis — there 
exist one more constant of motion given by ^U] 

with b and c constants coming from the parametrization, in cylindrical coordinates, of 
such surfaces as 

z = bp + c, (9) 

where / = Iuj is the angular momentum of the sphere respect to the normal to the 
surface of revolution at the point of contact and u the rotational angular speed around 
this same axis. In the conical surface b is the slope of the generatrix and c is the z- 
coordinate of its vertex. The easiest way of showing that (jSJ) is a constant is by taking 
its t-derivative and using to show that it vanishes. In is also worth pinpointing 
that z > > z = —a esc a since the sphere is geometrically constrained — i. e. we are 
assuming that it is not possible for the sphere to cut through the surface of the cone; z 
being the ^-coordinate of the vertex of the cone on which the sphere is rolling. 
Using equation (j2J and the expression for p in the integrals of motion (jHJ), fZJ) and (JHJ), 
we obtain 

1-15 5 

-M(f 2 + r 2 d 2 sin 2 a) + -Ma 2 u 2 + -Mgr cos a = -E, (10) 
2 T T T 

Mr 2 d sin 2 a = L z , (11) 

L 2 cosa 2L+ _ 

auj t, — = = n, (12) 

Mr sin 2 a 5 Ma 

where h is a new constant. Notice that equation (jllj) shows that never changes sign, 
i. e. the change in d is monotonic. 

4. The equivalent radial problem 

Combining the previous equations, we can get a new expression for r which can be 
expressed as 

£ = lMr 2 + 4 + - + CV, (13) 
2 r l r 

with £ , A, B and C constants defined by 

£ =^ E --Mh 2 , 
7 7 



L 2 (l + | cot 2 a) 
— - i 

2M sin 2 a 

n 2 L z h 

B = cot a, 

7 sin a 

5 

C = -Mucosa, 



(14) 
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Figure 2. The effective potentiai V„a{r) is plotted against r in arbitrary units. A 
particular value of the constant 'energy' £ = — 1 is also shown. The points marked 
n and ri are turning points, that is, points in which r = 0. As the plot exhibits the 
effective potential has only one minimum and thus just two turning points for any 
£ > V" c g-(r c ). The plot also exhibits that V^fj(r) attains its minimum value at r c and 
that this is a stable minimum. The plot was made using A = 1, B = — 3.5, C = 2 in 
equation (15). 



in the case the sphere is rotating around the cone axis [L z ^ and < a < ir/2] both 
A and C are strictly positive constants. 

Expression (fT3|) can be regarded as the energy equation for an equivalent, purely 
'radial' problem, with an effective potential energy term given by 

V eS (r) = £ + j + Cr, (15) 

which we show in figure 2. In this figure we have also traced a constant £(= 
—1, in arbitrary units) value as the horizontal line crossing V^g in two points, r\ and 
r2 — the turning points where r vanishes. This figure shows the typical behaviour of V e Q. 
A glance to figure 2 thus shows that the motion is restricted between the two turning 
points, meaning that the sphere's CM moves on a strip bounded by two parallel circles 
on the conical surface. 

The graph of V g q, figure 2, also shows a minimum value at r c . That is, an stable 
circular orbit is possible at this distance r c . The actual radius of the circular orbit, as 
measured from the cone's symmetry axis, is p c = r c sina. 

Moreover, equation (fTo^ can be directly integrated to yield the general solution of 
the problem as 

f r dr' 

t = ±l , (16) 



n 
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where we have chosen the reference point as one of the turning points of the motion. 
Taking into account that r = (dr/dt)d, we obtain, from equations (jllj) and ([13)1 . the 
orbit equation as 

* - ± ,V f , L ' d/ (it) 

sin «V2M J ri r / 2 /g _ (_A + s + Cr ,^ 

where r\ corresponds to the first turning point and the double sign corresponds to the 
clockwise and counterclockwise rotations, respectively. Equations (fTH|) and fP7|) are the 
general form of the solution of the problem in terms of quadratures. The integrals (|17|) 
and (jlfij) can be expressed in terms of elliptic integrals of the first and of the second 
kind, but such expressions do not offer much insight into the properties of the motion. 
However, a simple analysis of figure 2 says that the motion is always bounded between 
ri and r 2 , and that, for given E, h, and L z , there can be only one circular orbit which 
besides is stable. To determine the radius, r c , of this orbit we must take u = — which 
means that the acceleration, and thus the force, in the r-direction vanishes — in (jHJ) to 
get 



7$„rv sin a — 



5g — 2auiJ} c 



cos a = 0, 



the suffix c meaning magnitudes evaluated on the parallel circle r = r c . Since the first 
term in equation (fT£|) is always positive and cos a > because < a < ir/2, the circular 
orbit is only possible when the initial conditions are such that 

§ c u c < — . (19) 
2a y ' 

Equation (|18|) is quite useful since it allows evaluating any one of the three quantities 
r c , $ c , or uj c as a function of initial conditions of the remaining two for starting the 
motion on the circular orbit. If we know the initial values ujq and we could ascertain 
whether equation (fTHj) has a solution or not, just by checking if those initial values satisfy 
equation (|19|). Clearly equation (fTT?)) holds if u and $0 have contrary signs independent 
of their magnitudes. Now equation (|18|1 has one solution for uo c with initial conditions 
r and $0 while it has two solutions with initial conditions r and uj one of opposite 
sign to ujq and another of the same sign, but satisfying equation (|19)1. 
Equation |T8|) in terms of the constants of motion L z and h is 

hrl M 2 gsm 2 acosa-2r c ML z hsm 2 acosa-7L 2 z ^1 - j cos 2 a^j =0.(20) 

This is a cubic equation for r c that has always a real positive root; but let us pinpoint 
that, due to condition (fT9*j) . this equation is only valid when 2hL z < 5Mr 2 gsin 2 a. The 
existence of only one real and positive root, follows from the fact that the last term in 
(j20j) is less than zero, whereas the coefficient of the first term is positive jT3] . However, 
it can be much simpler to just take a look at figure 2, as it exhibits that there is one 
and only one circular orbit on the surface. 



Sphere rolling on the surface of a cone 



8 



There is a further point worth mentioning, the existence of a relationship between 
the radius of the circular orbit and the r- values of the turning points, r i5 i = 1,2, 
namely 

, 2r 2 r 2 B ( rir 2 \ , „ 

This is the sort of property called universal in in fact, (|21|) is very easily shown to 
reduce to equation (13) of reference jB] in the point particle limit (that is, when B = 0). 

5. Analogy with an Non Homogeneous Harmonic Oscillator and the Orbit 
on a Preelly Falling Cone. 

Let us substitute equations (0) in the first of equations (JSJ), to get 

2 <p . . 5 

r — rep — (aoM cos a — -gcosa, (22) 

7 sin a 7 

where ip = i? sin a. Multiplying by M and using equations© and (JSJ) in the right hand 
side of equation (|22[). we get 



where 



and 



M(f-r0 2 ) = /(r), (23) 



/(r) = ^ + ^-^oosa, (24) 



/ 2 
= -Ln cot a 
7 

,. 2 L 2 cot 2 a ,„ r > 

A = (25) 

7 M 

2, ' 



L = Mr l (p = 

sin a 

therefore the equations of motion of the sphere are found, as we did in section 4, 
analogous to the equations of motion of a particle in a central force field. The variables 
r(t) and <p(t) can be considered as the polar coordinates of a particle moving on a plane 
under the action of the "central" force f(r) given in equation (|23j). If, additionally we 
assume the cone is in free fall, that is, that g = 0, then the orbit of the center of the 
sphere, r = r (if), can be easily obtained by rewriting equation ()22j) in terms of ip, using 
the relations 

d . d L d 

Jt =(f d^ = Mr^d^ ^ > 



and 

d 2 L 2 d ( 1 d 



dt 2 Mr 2 dip \r 2 dip 



(27) 
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Now, introducing Binet's variable W = 1/r, the equation of motion (J23j) becomes, with 
# = and using (JH| and (j2Zj), 

+ l + -cot 2 a )W=-- (28) 



d(p 2 V 7 / 7 L 

or, in terms of the angle in the polar plane, $ = sin a, and the vertical component of 
the angular momentum, L z , 

d 2 W / 5 9 \ TTr 2 M/i cos a sin 2 a , . 

+ I 1- -cos 2 a ) VP= . (29) 



dtf 2 V 7 / 7 L z 

This is the well-known inhomogeneous differential equation for harmonic motion. If we 
define fl 2 = (1 — (5/7) cos 2 a), the general solution can be written 

1 . „ 2Mh sin 2 a cos a , , 

W = - = ^Icos - O )] ^5 • (30) 

The constants of integration A and ?9o can be determined by the initial conditions. 
To clasify the solutions of the problem in free fall, we recognize two cases, 

A) when h and L z have different signs, i.e. hL z < 0; or 

B) when h and L z have the same sign, i.e., hL z > 0. 

In any case, the general solution can be written in the form 

I ens rv 

W = - = — — ecos[ft(tf-tf )±l] (31) 
r V 



with 

7\L z \n 2 

2M\h\ sin 2 a 

or in the form 



P = W77iTT—2Z ( 32 ) 



with 



e ± 1 

r(&) = r min — — — (33) 

e cos \l(v — v ) ± 1 



(34) 



cos cue ± 1 



where e and $0 are constants of integration. The sign in equation is plus if hL z < 
and minus if hL z > 0. Without loss of generality we can assume that e > 0, since its 
sign can be changed by substituting •& — > d + 7r, that is, just by changing the orientation 
of the coordinates. 

The possible motions are best analysed case by follows. 

A) In the case that hL z > 0, we have the minus sign in equation ()34j) and therefore 
e > 1. In this case there are two values of •& for which e cos f2(i9 — $0) — 1- When d 
approaches any of these two values r — > 00; therefore they define two asymptotes 
and the orbit corresponds to a branch of an hyperbola. 
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B) In the case that hL z < 0, the plus sign is in order in equation (|34j) . We find now 
four possibilities, e = 0,0<e<l,e = l and e > 1. The motion in these cases is 
as we describe in what follows. 

Bl. If e = 0, the orbit is the circle 

r(i?) = r c = p esc a — - cos 2 c^j ; (35) 

this relation coincides with (1201) if g = 0. 
B2. If < e < 1, the orbit is confined to the strip defined by the two parallel 
circles on the cone, having radii 

r = r min (36) 

and 

l + e 



mm -, 

1 — e 



(37) 



The CM of the sphere touches the circle with radius r = r max when 

$-$ = ^ = 71 (38) 
| cos 2 a 

which means that the apsidal angle ip (the angle swept by the CM's radius 
vector in going from r min to r max ) is 

^ = ! > 7T, (39) 

1 — | cos 2 a 

thus the apsides advance in each period of rotation. As the orbit is symmetric 
with respect to the polar axis (r(— i?) = T (■$)), the angular displacement of the 
apsides in each rotation of the center of the sphere around the z-axis is 

Aip = 2tt[(1 - j cos 2 a)' 1 ' 2 - 1] (40) 

and the angular velocity of precession of the apsides is 

^ = ^[(l-^cos 2 ^.!] = W[(l-^cos 2 a)- 1 / 2 _l](41) 

where T is the period of rotation of the CM around the z axis, and 
1 f T A , 2tt 



(0) = ^ J Q Mt=— (42) 
is the mean value of 

Since (1 — 5 cos 2 a/7) 1 / 2 can be shown to be an irrational number for < a < 
tt/2, the orbit never repeats itself thus filling completely the strip that contains 
the orbit as t — > oo. We may say that the orbit is ergodic on the strip, that is, 
that eventually the orbit would densely fill all of the strip's surface. 
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and therefore r — > oo when d approaches any of these values, which simply 
define the asymptotes of the corresponding branch of the hyperbola. 



The only thing that remains to be done is to obtain the value of e in terms of the 
constants of motion. For this we have only to susbstitute the values of r min and of p 
into equation (JT3J) (with C = 0) to get 



Note that e depends only on £ and h, and, furthermore, that < e < 1, if £ < 0; that 
e = 1, if £ — 0; and that e > 1, if £ > 0. But, if hL z > then B > and therefore 
£ > 0, from which we have e > 1 always, consequently this case is equivalent to motion 
in a repulsive force field. 

6. Conclusions 

We have studied the motion of a sphere rolling on the inner surface of a right circular cone 
with opening angle a under the assumption that it does not slip. The motion was first 
considered under the action of gravity and we managed to obtain an implicit solution in 
quadratures. We established the existence of bounded motions and of circular orbits in 
the system analysing an effective potential in which the equivalent single particle system 
moves. With the help of this effective potential we were able to reduce the problem to 
quadratures obtaining implicitly r(t) and thus 6{t). 

However, to actually integrate the equations of motion for obtaining the centre's 
of mass orbit in terms of elementary functions, we required changing to a freely falling 
frame. In such a frame we have been able to obtain solutions for the CM orbit analogous 
to ones in a related central problem. We have been able to show, with the help of this 
relation, that the CM can move on hyperbolic or on preceding elliptic, but never on 
parabolic orbits. 

The analogy between motions in our system with central-field motion and specially 
with an inhomogeneous oscillator in the freely falling case, should be pointed out as 
interesting and quite amusing. This aspect of the behaviour has a certain pedagogical 
significance as can serve to emphasise that the same equations have always the same 
solutions; it does not really matter that we are describing different physical systems. 




(46) 
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